Mult i- quantum eigenstates of a linear chain of coupled qubits 
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We present a technique to identify exact analytic expressions for the mult i- quantum eigenstates 
of a linear chain of coupled qubits. A choice of Hilbert subspaces is described which allows an exact 
solution of the stationary Schrodinger equation without imposing periodic boundary conditions and 
without neglecting end effects, fully including the dipole-dipole nearest-neighbor interaction between 
the atoms. The treatment is valid for an arbitrary coherent excitation in the atomic system, any 
number of atoms, any size of the chain relative to the resonant wavelength and arbitrary initial 
conditions of the atomic system. The procedure we develop is general enough to be adopted for 
the study of excitation in an arbitrary array of atoms including spin chains and one-dimensional 
Bose-Einstein condensates. 
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It is well established that a large scale quantum com- 
putation will require a large number of strongly coupled 
atoms (qubits) and in particular multi-photon excitations 
to create large quantum superpositions of different quan- 
tum states of the system flj . Naturally, this means that 
an operation of such systems will require collective exci- 
tations involving more than one photon. Although collec- 
tive excitations in multi-atom systems have been studied 
in the past, it was typically done in the context of ra- 
diation properties of a large number of atoms randomly 
distributed in space, such as in a gas cell 0- There have 
been many exact studies of collective eigenstates of spa- 
tially ordered systems, such as Heisenberg chains 0,0>l|- 
However, these studies suffer from one common draw- 
back: they are limited to finding the ground state or 
one-photon excited states in the limit of a large number 
of atoms, or are subject to periodic boundary conditions. 

With the recent progress in trapping and cooling of a 
small number of atoms or ions, the attention has been 
drawn towards systems that are comprised of a small 
number of atoms having a definite geometrical arrange- 
ment, as dictated by the confining field of a hypotheti- 
cal atom trap 6]. A number of analytical studies have 
been performed on systems composed of N atoms con- 
fined to fixed positions, with a particular interest in ring 
arrangements [?J, perhaps as they have periodic bound- 
ary conditions which makes the eigenstates easier to ob- 
tain. However, exact analytic studies, which give ex- 
plicit forms of the eigenstates, have been performed on 
systems containing only two and three atoms 0, H, @ • 
More gen eral results valid for arbitrary N have also been 
given |ld , but are limited to one-photon excitations 
only . More exotic shapes, such as diamond structures 
|12| . have also been investigated, but the calculations 
were limited to two-photon excitations and small num- 
bers of atoms. While these schemes successfully demon- 
strate the situation for obtaining one- or two-photon- 
excited states in A^-atom systems, the analysis for the 
excitation of such a system by an arbitrary number of 
photons, giving a multi-quantum eigenstate, is unknown. 



Other treatments of these systems stem from statistical 
physics, where properties are extracted from the partition 
function, which is determined by the eigenvalues of the 
system. The partition function can provide a witness to 
entanglement |13|. it does not, however, give information 
about the energy eigenstates of the system. These states 
may have applications in quantum information and need 
to be quantified in order to perform quantum computa- 
tion. This forms the motivation of this work to derive 
the explicit analytic expressions for the multi-quantum 
eigenstates of a linear chain of A" atoms that can be ap- 
plied to a variety of settings. The availability of these 
eigenstates is highly advantageous because it provides a 
convenient ground for rigorous examination of entangled 
properties of a particular arrangement of atoms and fa- 
cilitates the explicit study of the stability of entangled 
states, something which can be obscured if we remove 
this generality. It may also provide an interesting tool 
to study entanglement creation in ID Bose-Einstein con- 
densates ^4|. We cover the case of the linear chain as 
although it is one of the simplest structures that can be 
constructed in an atom trap, it has not been solved ex- 
actly while imposing the generality that we demand here. 
There is also a considerable interest in analyzing multi- 
atom systems to realize a phase gate operation [15J, in 
which we do not impose periodic boundary conditions 
and include the end effects. 

In this paper, we study the creation of multi-quantum 
eigenstates in a linear chain of A" identical, equally spaced 
and confined to fixed positions, two-level atoms each in- 
teracting only with its nearest neighbors through the 
dipole-dipole interaction. This model might be realized 
in practice by placing the atoms in a tightly confi ning lin- 
ear trap, an optical lattice or in an atomic chip |161 Il7|. 
The Hamiltonian for the linear chain of atoms is given by 
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h = H Q + v = huj J2s* + nn ^s~, (i) 
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where Hq is the interaction-free Hamiltonian and V is 
the dipole-dipole interaction between the atoms. 

In Eq. 0}, w o is the transition frequency of a two- level 
atom in isolation, Sf is the energy operator of the i-th 
atom, and Sf~ , S~ are the raising and lowering operators 
for the i-th atom, respectively. The dipole-dipole inter- 
action parameter f2 = O(ry-) depends on the distance 
between adjacent atoms in the linear chain and the po- 
larization of the atomic dipole moments relatively to the 
interatomic axis 

To find the multi-quantum energy states of the system, 
we have to solve the stationary Schrodinger equation with 
the Hamiltonian H. We propose a procedure which, de- 
spite of the complexity of the problem, facilitates an ex- 
act analytic solution of the Schrodinger equation valid for 
an arbitrary number of excitations M. The interaction- 
free Hamiltonian Hq has N + 1 energy levels of energies 
E ( M) = Mhuo, where M = 0, 1, 2, . . . , N. Since there 
are N two-level atoms, it follows that Ho operates on 
a 2 Ar -dimensional Hilbert space W. We can index the 
atoms in the chain by the numbers 1 to N; collective ex- 
citations the system can then be represented using the 
ket \ki, ■ ■ ■ , kit), where the non-zero integers k\, . . . , &m 
denote the indices of the atoms which are in their excited 
state: 

\k 1 ,...,k M )=S+---S+ M \0), (2) 

where |0) denotes the ground state of the collective sys- 
tem, and is nondegenerate. The ground states |0) and 
excited states \k\, . . . , &m) form a basis which spans W. 
The maximally excited interaction-free energy eigenstate, 
i.e. the collective state where all N atoms are in their 
upper state, is then represented by |1, 2, . . . , N) in our 
notation. 

Under the action of Hq, a particular energy level M 
has a degeneracy factor of N\/(N — M)\M\, which is 
equal to the number of ways N atoms can be divided into 
two groups, with M atoms excited and N — M atoms in 
their ground states. However, while H acts on the same 
Hilbert space W, it leads to a significant splitting of the 
degeneracies of the collective system. 

We observe that the process of finding the eigenstates 
of H may be simplified by noting that Ho and V com- 
mute. Since 

\H V] N 

i-^i= J2 (Sf[S?,S^ + [S?,S+]S^, (3) 

k=l 

li-fe|=l 

and using the relations [Sf, S^] — ±S^Sij, we can easily 
see that [Ho, V] = 0. This result, coupled with the fact 
that Hq and V are Hermitian means that the operators 
Ho and V share a complete set of orthonormal eigenstates 
{!<&)}• Since H is just the sum of these two operators, 
and {|</>i}} is a complete set, it follows that H, Hq and V 
share the same complete set of orthonormal eigenstates 



{!</>»)}• The Ho operator allows us to partition the state 
vector space W into N + 1 subspaces W- 7 : 

N 

W=($W j , (4) 

1=0 

where the subspace corresponds to all states which 
give an eigenvalue of jfiui when operated on by Ho, and 
is of dimensionality Nl/(N — In this way, we see 

that Hq is an automorphism T(Hq) : W j — > for all 
W J . Since H and V commute with Ho, they must be 
mappings of the form — > W J for all subspaces WK 
Thus, to find the energy eigenstates of the system, we can 
solve the Schrodinger equation for H for each subspace 
Wi instead of W . However, all states in a particular 
are eigenstates of H , so we can neglect this operator 
and simply solve the Schrodinger equation for V. To 
accomplish this, we take the most general state vector 
one can form in the subspace W M : 

|^ (M) ) = C(k u ...,k M )\ki,...,k M ), (5) 

ki <---<kM 

where the summation runs through all M-tuples 
(ki, . . . ,k M ) such that 1 < k± < ■■■ < kni < M. To 
then obtain the energy eigenstates, we solve 

V\^) = AE\^), (6) 

where \ip^ M ^) S W M and AE is the energy level shift as- 
sociated with the dipole-dipole interaction V. The eigen- 
value associated with H is then given by 

H\^) = Ho\^ M) ) + V\^) 

= {E ( M) +A J E)|V (M) ). (7) 

We assume that the energy eigenstates of the inter- 
action Hamiltonian are of the form JSJ and constitute 
a complete set. In order to determine the coefficients 
C(fci, . . . , ku), w e substitute 10 into Eq. @J, and 
find the following multi-term recurrence relation relating 
the coefficients C(k\, . . . , /cm) to the energy eigenvalue 
shift AE: 

M 

AEC(k 1 ,...,k M ) = mY^[C(ki,...,k j + l,...,k M ) 

3=1 

+C{k 1 ,...,k J -l,...,k M )}, (8) 

where it is assumed that the summation in the above 
expression must not include terms which have an invalid 
value for kj, i.e. are not generated by Eq. JfjJ). There are 
only three such cases: (a) when k\ is equal to zero or (b) 
ku = N +1 and thus the index cannot refer to an actual 
atom, and (c) when two of the kj are equal, which would 
rule out the state from being an eigenvector of V. 

Equation © without such a restriction on summa- 
tion can be solved easily, so we prescribe the following 
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method: we allow invalid terms to be included in Eq. 
(|8*|) , on the condition that they must vanish. This means 
that we now regard cases (a), (b), and (c) as conditions 
under which a coefficient must vanish. 

By inspection, we see that Eq. © resembles a pattern 
that is seen in sine or cosine functions f(x), for /(x + l) + 
f(x — l) oc f(x). As a trial solution, we thus consider the 
general form 

C(ki, . . . , k M ) = K a ... d j a {a a k 1 + A a ) ■ ■ ■ 

■■■f d (S d k M + D d ), (9) 

where a a , . . . , 5 d and A a , . . . , D d are arbitrary constants 
along with K a ... dl and the summation range for the M 
indices a, . . . ,d ranges from 1 to M. By operating on the 
kj individually by sine or cosine functions, this expansion 
can be made fully consistent with Eq. (jSJ) using further 
restrictions which now follow. 

The only way that condition (c) can be satisfied is if 
K a ... d is antisymmetric in any two pairs of indices. Con- 
dition (a) then implies that cither all the fi are sine func- 
tions with the offsets Ai, . . . ,D d equal to zero or the fi 
are cosine functions where Ai,...,Dd are all equal to 
7r/2. These solutions are not independent - they are 
equal up to a multiplicative constant, so we choose the 
sine function. The antisymmetry property also requires 
that oti = ■ ■ ■ = Si. We then have 

C(fci, . . . ,k M ) = Ke a ... d sm{a a ki) • ■ ■ Bin(a d k M ), (10) 

where e a ... d is the permutation (Levi-Civita) symbol and 
summation is implied over the repeated M symbols 
a, . . . ,d from 1 to M. The constant K will be set equal 
to unity from now on as it is not constrained by any 
conditions. 

For condition (b) to be satisfied, and owing to the per- 
mutations of the di, the constants a\, . . . , cxm must equal 
g\9, . . . , <7a/0, where 9 — n/(N + 1) and the gi are inte- 
gers. Our solution now becomes 

C(fci, . . . , k M ) = e a ... d sva(g a ki9) ■ ■ ■ sm(g d kM9), (H) 

or equivalently 

C(fci, . . . , k M ) = e a ... d sm(gik a 9) ■ ■ ■ sm(g M k d 9). (12) 

We have thus derived the coefficients apart from the 
arbitrary integers gi,..., gu- Clearly, no two gi can be 
equal, for the C(k%, . . . , fc^r) would vanish due to the 
presence of the permutation symbol. Also, to avoid du- 
plicate cases, we require for non-trivial eigenstates that 

9i < ■ ■ ■ < 9m- (13) 

If a particular gi is equal to N + 1, the coefficient 
C(fci, . . . , fcjv/) vanishes; therefore all the gi must be less 
than or greater than N + l. The difference between these 
two cases amounts to multiplying Eq. (|12|) by a minus 



sign on the right-hand side. Since such a factor is irrel- 
evant, we impose the restriction that gi < N for all gi. 
Due to the properties of the sine function, this implies 
also that gi > 1 for all gi. Thus we can state constraints 
on the gi\ 

1 < 31 < • • • < 9m < N. (14) 

Now that we have determined the coefficients of the 
eigenstates lf5fk we are in a position to determine the 
energy eigenvalue shift. Using JSJ, we have 

AEC(kt,...,kM) 

M 

= hfl^22 cos{g c 6) [e a ... c ... e sin(g a ki9) ■ ■ ■ 

3=1 

• • • sin(g c kj9) ■ ■ ■ sin(g e k M 9)} , (15) 

where we have converted pairs of terms of the form 
sing c (fcj + 1)9 and sing c (fcj — 1)9 into the product 
2 cos g c 9 sin g c kj9. We repeat this for each sine pair, giv- 
ing the end result 

^EC(k u ...,k M ) = 2nnc(k 1 ,...,k M ) 

M 

x^cosM). (16) 

i=i 

The energy eigenvalue shift is thus given by 

M 

AE = 27ifi^cos(pi0). (17) 

i=l 

We now collect the results to give the energy eigen- 
states and eigenvalues for the system Q, valid for an 
arbitrary number of atoms N and an arbitrary number 
of excitations M. The unnormalized energy eigenstates 
of (JIJ are given by 

\4^-9m)= E C^lh,...,^), (18) 

fel<-"<&M 

with corresponding eigenvalues 

M 

e^ ] . bu = Mnw + 2nnJ2 cos (9iO) (19) 

i=l 

where 

Cgl-'-gM = M9aki9) ■ ■ ■ sm(g d k M 9), (20) 

and 9 = ir/(N +1). The parameter M equals the num- 
ber of excited atoms in a particular eigenstate, and thus 
forms one of the quantum numbers that enumerate the 
energy eigenstates. The other quantum numbers are the 
M numbers gi, which are subject to the constraint 114fl . 
There are N\/{N — M)\ ways of choosing values for the gi 
so that for all g%, 1 < gi < N. However, the con- 
straint (|14|) forces us to divide the number of combina- 
tions by the number of ways one can re-label the gi, which 
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is M\. Hence there are N\/(N — M)\M\ energy eigen- 
states associated with a given value of M. Note that the 
number of the states for a given M is the same as it was 
in the interaction- free system mentioned earlier. 

To fit the ground eigenstate (M = 0) into the solu- 
tion Q18JI. we observe that C 1 ^ 1 .'.'.'^ has no gi's and no fcj's 
for the case M = 0, so we define C* 1 = C, an arbi- 
trary constant which may be set equal to unity. Clearly 
the summation sign in (|18|l disappears in this case, and 
the only ket present is defined to be that corresponding 
to the ground state of the system, |0). 

One can see the step which enabled us to determine the 
exact analytic expressions for the multi-quantum eigen- 
states. The recurrence relation derived from the eigen- 
value equation had restrictions on its summation range; 
however, it is clear that there is a simple pattern in the 
summation formula. We completed the pattern by in- 
cluding invalid terms, and created boundary conditions 
which made these extra terms vanish. Other treatments 
handle the problem of the restrictions on the original 
summation by imposing periodic boundary conditions. 
However, the latter treatment implicitly converts the lin- 
ear chain into a ring - a different physical system. That 
the derived states are eigenstates of a linear chain Hamil- 
tonian is the justification of our approach. 

As an illustration of our general solution, we discuss 
two specific examples, M = 1 and M = 2, corresponding 
to single- and two-photon excitation of the iV-atom sys- 
tem. The general eigenstates <|18[) for M = 1 simplify to 

N 

\^) =J2M9ik6)\k), (21) 

k=l 

where 1 < g\ < N is a single quantum number differen- 
tiating between N states. The associated eigenvalues are 

E$ = Hlu + 2MI cos gi 6. (22) 



For these states, we can see an interesting property in the 
small sample model, where the interatomic sep aration is 
much smaller than transition wavelength [l||. For the 
states where g± is even, the transition dipole moment 
{ipg 1 |/i|0), where |0) is the ground state, vanishes. The 
reason for this is that the coefficient of \k) is the negative 
of the coefficient of \N + 1 — k) in the above eigenstate 
expansion for such states. The sum of these coefficients 
then gives zero for the transition dipole moment. Thus 
gi-even states do not spontaneously decay to the ground 
state. 

For M = 2, the eigenstates are of the form 

JV-l N 

A) = E E [sin(ffiM)sinG&M) 

ki — 1 fe2=fei+l 

- sin(g 1 k 2 9) sm(g 2 k 1 9)]\k 1 , k 2 ), (23) 

with energies 

E(f g2 = 2Hlo q + 27LO(cos g x 6 + cos g 2 6). (24) 
For this case, 1 < gi < g 2 < N, giving a total of 
N(N — l)/2 eigenstates. This is the first excitation level 
where the antisymmetric property of the Cg^.-.g*/ ap- 
pears. For a small number of atoms, N = 3, the above 
eigenstate and eigenvalue examples agree (apart from the 
normalization constant) with the case of the isosceles tri- 
angle flattened into a straight line, with the dipolc- 
dipole interaction between the two atoms at the ends of 
the line removed. 

In summary, we have proposed a simple and ro- 
bust method to derive exact expressions for the multi- 
quantum states of a linear chain of N trapped atoms 
interacting with its nearest neighbors through the dipole- 
dipole interaction. The method is based on a technique 
of converting the eigenvalue equation into a specific re- 
currence relation, rather than the usual technique of di- 
agonalizing a matrix, the latter being more suited to the 
case where the number of atoms is not arbitrary. 
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